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Abstract 

Searches for hght charged Higgs bosons {H ) in the decay of top quarks, t — )■ H b, are being 

r^ ', carried out at the LHC and at the Tevatron. It is assumed that the dominant decay channels for 

II I such an H^ state are either H^ — )• tu or H^ — )• as, and separate searches are performed with 

^H' comparable sensitivity to the parameters mjj± and tan /3 of the scalar potential. The branching 

r~| ' ratio for the decay H^ — )• cb can be as large as 80% in the Aligned Two Higgs Doublet Model and 

in models with three or more Higgs doublets with natural flavour conservation, while satisfying 

CN ! the constraint from 6 — )■ 57 for m^/i < rrit- Although the current search strategy for H^ — )• as 



is also sensitive to H — )■ cb, a considerable gain in sensitivity could be obtained by tagging the 



VO ' b quark from the decay H — )■ cb. Such an analysis, which could be readily performed at the 

Tevatron and in the 7 TeV and 8 TeV runs of the LHC, would probe a parameter space of the 
^^..^ fermionic couplings of H in the above models which at present cannot be probed by experimental 

f^ , observables in flavour physics. 
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I. INTRODUCTION 

There is much ongoing experimental effort by the ATLAS and CMS collaborations at 
the CERN Large Hadron Collider (LHC) to search for the neutral Higgs boson (h^) of the 
Standard Model (SM) [l|, l2|. This model of spontaneous symmetry breaking will be tested 
at the LHC over all of the theoretically preferred mass range, in an experimental programme 
which is expected to be completed by the end of the 8 TeV run of the LHC. At present (Sj 
there are only two regions for the mass for h^ which have not been excluded at 95% c.l: i) 
a region of light mass, with 122 GeV < rriho < 128 GeV, and ii) a region of heavy mass, 
niho > 600 GeV. 

However, this simplest model of one fundamental scalar with a vacuum expectation value 
(vev) might not be nature's choice. There could be additional scalar fields which also 
contribute to the masses of the fermions and weak bosons, with a more complicated scalar 
potential which depends on several arbitrary parameters. Importantly, even in the event 
of no signal for a SM-like Higgs boson at the LHC, the search for scalar particles should 
continue in earnest due to the fact that a non-minimal Higgs sector can give rise to different 
experimental signatures, some of which are challenging to detect. Consequently, it will take 
much longer for the LHC to probe all of the parameter space of such models. 

A commonly-studied extension of the Higgs sector of the SM is the "Two Higgs Doublet 
Model" (2HDM), which is composed of two Higgs isospin doublets J4| (this model has recently 
been reviewed in |5|). Notably, this structure is necessary in the Minimal Supersymmetric 
(SUSY) extension of the SM (called the "MSSM"). The extra Higgs doublet gives rise to a 
particle spectrum with multiple Higgs bosons; three are electrically neutral (two are CP-even, 
one is CP-odd) and two are electrically charged (denoted by H^). Flavour-changing neutral 
currents (FCNCs) mediated by scalars at tree level can be eliminated by requiring that 
the scalar interactions with the fermions are invariant under discrete symmetries ("natural 
flavour conservation" , NFC) |6| . The discovery of a charged scalar H"^ would be unequivocal 
evidence of a non-minimal Higgs sector, and there have been many studies of the prospects 
of directly observing if^ from a 2HDM or the MSSM at the Tevatron and LHC [7(] (for 
reviews see p, l8|) Moreover, the effect of H^ on the decay rates of mesons (especially B 
mesons) also plays a major role in constraining r/i/fi and the fermionic couplings of H"^ 

0.0. ^ 

If mH± < mt + nib, such particles would most co pio usly (though not exclusively [lO(]) be 
produced in the decays of top quarks via t — )■ H'^b [ll|. Searches in this channel have been 
performed by the Tevatron experiments, assuming the decay modes H^ — )■ cs and H^ — )■ tu 



|12l . |13| . Since no signal has been observed, constraints are obtained on the parameter space 
of [m//±, tan/3], where tan/3 = ^2/^1 (i.e. the ratio of the vacuum expectation values of the 
two scalar doublets). Searches in these channels have now been carried out at the LHC: i) 



for H^ -> cs with 0.035 fb"^ bvATLAS [ij, and ii) for H^ -^ tu with 4.8 fb'^ by ATLAS 



15[ and with 1 fb~^ by CMS [16|. These are the first searches for H"^ at this collider. The 
constraints on [mH±,ta.nP] from the LHC searches for t — )■ H'^b are now superior to those 
obtained from the corresponding Tevatron searches. 

The phenomenology of if ^ in models with three or more Higgs doublets (called Multi- 
Higgs Doublet Models, MHDM), was first studied comprehensively in [17|, with an emphasis 
on the constraints from low-energy processes (e.g. the decays of mesons). Although the 
phenomenology of H"^ at high-energy colliders in a MHDM and in a 2HDM has many 
similarities, the possibility of mH± < fnt together with an enhanced branching ratio (BR) 



for H"^ — )■ ch would be a distinctive feature of the MHDM. This scenario 
mentioned in 17|] and studied in more detail in |l8l420l 
ongoing searches for t — ?■ H^^h with H"^ — )■ cs by the LHC [14 



which was first 

is of immediate interest for the 

Although the current limits 



on H^ — )• cs can be applied to the decay H^ — >■ ch (as discussed in [2]| in the context of the 
Tevatron searches), a further improvement in sensitivity to t — >■ H^h with H^ — > ch could 



We will estimate 
ch) and to the fermionic couplings of H"^ in this 



be obtained by tagging the h quark which originates from H^ 

the increase in sensitivity to BR(i7^ 

scenario. 

Large values of BR(if ^ — )■ ch) are also possible in certain 2HDMs, such as the "flipped 
2HDM" with NFC [l8|, |2l|, |23| . However, in this model one would expect mH± > rrit due to 
the constraint from h ^ sj {mH± > 295 GeV |24|-|26l|). and thus t — )■ H'^h with H^ — )■ ch 
would not proceed unless there were additional New Physics beyond that of the 2HDM which 
contributed to 6 — )■ S7, and weakened the constraint on mH±- In the "Aligned Two Higgs 
Doublet Model" (A2HDM) Q there are no FCNCs (as is the case in a 2HDM with NFC) 
due to an alignment of the Yukawa couplings. The phenomenology of H^ in the A2HDM 
[28| is very similar to that of if ^ in a MHDM 



30 



and mjj± < rrtt in the A2HDM is also 
compatible with constraints from h — )■ 57. Our numerical results for BR(ilf^ — >■ ch) in a 
MHDM apply directly to the A2HDM. The 2HDM without NFC and without alignment 
also has a sizeable parameter space for a large 'BK{H^ — )■ ch), and a detailed study can be 
found in |22l ]. 

Our work is organised as follows. In section II we describe the fermionic interactions 
of H"^ in the MHDM/A2HDM. In section HI we quantify the parameter space for a large 
BK{H^ -^ ch) in the MHDM/A2HDM. In section IV we summarise the Tevatron/LHC 
searches for t — )■ H^h with H^ -^ cs and discuss how they could be optimised for H^ — )■ ch. 
Section V contains our numerical results, with conclusions in section VI. 



II. H^ IN MODELS WITH MORE THAN TWO HIGGS DOUBLETS AND IN 
THE ALIGNED 2HDM 



In a general 2HDM each fermion type (i.e. up-type quarks, down-type quarks and charged 
leptons) couples to both of the scalar doublets. This would lead to FCNCs which are 
mediated by the neutral scalars, and the magnitude of the associated Yukawa coupling is 
constrained by experimental data (especially meson-antimeson mixing and the decays of 
mesons). Such FCNCs can be suppressed by assuming that the flavour-changing Yukawa 
couplings are very small, which can be achieved by invoking a specific structure of the fermion 



mass matrices 31 



An alternative approach is to eliminate the FCNCs by requiring that 
the Lagrangian is invariant under a discrete symmetry, which is achieved if each species of 
fermion couples to at most one scalar doublet (NFC). This condition leads to four distinct 
types of 2HDMs which differ in their Yukawa couplings. These four models are called 
Model I, Model II, Lepton-specific and Flipped.^ Models I and II have received much 
phenomenological attenti on, while the study of the Lepton-specific and Flipped models has 
been revived recently |2l|, |23|, |33|, |3J| , with early studies in [l8|, |20|, |35| . We now introduce the 
fermionic couplings of H^ in a MHDM and the A2HDM, and discuss how these couplings 



^ The Legton-spccific and Flipped models are referred to as IIA and IIB in [S^l , IV and III in [7| , I' and II' 



in 17-20i and X and Y in 



23 



differ from those in the above 2HDMs with NFC. 



A. The Multi-Higgs Doublet Model (MHDM) 



A MHDM is an extension of the 2HDM with n scalar doublets, where n > 3. The 
suppression of FCNCs is obtained by imposing NFC. As in the 2HDM, the MHDM has the 
virtue of predicting p = 1 at tree level, with finite higher-order corrections which depend 
on the mass splittings of the scalars. In the MHDM there are n — 1 charged scalars, and 
a detailed study of the phenomenology of the lightest if ^ in such models was performed 



m 



17 



with the assumption that the other H"^ are much heavier. The interaction of the 
lightest H^ in a MHDM with the fermions is described by the following Lagrangian: 



£ 



v^K 



H± 



ud- 



-U 



[mdXPR + muYPL)dH' 



y/2me 



Z^l^rH^ 



H.c. 



(1) 



Here u and d denote up-type quarks and down-type quarks respectively (for all three 
generations); Vud is a CKM matrix element; m„, md and nif, are the masses of the up- 
type quarks, down-type quarks and charged leptons respectively; Pl and Pr are chirality 
projection operators, and v = 246 GeV. In a 2HDM with natural flavour conservation, the 
couplings X, Y and Z are determined solely by tan/3 = f2/fi. The values of X, Y and Z 
in the four versions of the 2HDM J7| are given in Table [H It is clear that \X\, \Y\ and \Z\ 
are simply related in the 2HDM e.g. one has \X\ = \Z\ = l/\Y\ for the Type II structure. 





X 


Y 


Z 


Type I 


-cot/3 


cot/3 


-cot/3 


Type II 


tan/3 


cot/3 


tan/3 


Lepton-specific 


-cot/3 


cot/3 


tan/3 


Flipped 


tan/3 


cot/3 


-cot/3 



TABLE I: The couplings X,Y and Z in the Yukawa interactions of H in the four versions of the 
2HDM with natural flavour conservation. 



In a MHDM the couplings X, Y and Z are arbitrary complex numbers, which are defined 
in terms of the n x n matrix U which diagonalises the mass matrix of the charged scalars: 



X- 



U, 



di 



u. 



Yi 



U,. 



dl 



Uul 



z. 






(2) 



We follow the notation of 17| in which i = 1 corresponds to the couplings of the charged 
Goldstone boson, and i runs from 2 to n for the physical charged scalars. The fermionic 
couplings of the lightest H^ in a MHDM are taken to be X2, I2 and Z2. The subscripts d, u, e 
take any integer value up to n, and specify which of the n doublets couples to which fermion 
type e.g. for a Type II structure one sets d = e = 1 and u = 2, while for the "democratic" 
3HDM (e.g. |30|) one has ci = 1, u = 2 and e = 3. In a 2HDM, f7 is a 2 x 2 matrix with 



elements given by sin/3 and cos/3 (i.e. one free parameter). In a 3HDM, f/ is a 3 x 3 matrix 
with four free parameters which can be taken as tan/3 = Vujvd-, tan 7 = \l {v\ + v^jvc-, a 
mixing angle Q for the two B.^^ and a complex phase b. An explicit form of the matrix f/ 



for the 3HDM is given in 30 1. 



Due to the unitarity of the matrix U one can derive the following identities [17 



Y,X,Y: = l{ioYd + u), (3) 

i=2 
n 

Y,X,Z: = -l (forrf^e), (4) 

=2 

n 

^F,Z; = l(forw^e), (5) 



i=2 



and 



on on o 



$:iA-.p=^-i.$:iy.i^=^-i,$:iz.i^=^-i, (6) 



".i ..2 



In a 2HDM these identities reduce to simple trigonometric relations involving tan/3. It 
is evident that the branching ratios of H^ to fermions in the MHDM depend on the three- 
dimensional parameter space of Xj, Yi and Zj, in contrast to the case in the 2HDM where 
a single parameter (tan/3) determines these three couplings. It is conventional to consider 
the phenomenology of the lightest if ^, assuming that the other if ^ are heavier. One then 
drops the i subscript on the couplings of the lightest if ^ and uses the Lagrangian in eq. ([T]). 

Many experimental observables in flavour physics would receive a contribution from if ^, 
and thus the magnitudes of X, Y and Z are constrained. Detailed studies have been per- 
formed in [l7| (in the context of a MHDM) and more recently in [28| in the context of the 
A2HDM (see also |9|). For mH± < mt these constraints are roughly as follows: |y| < 1 
from Z — 7- 66 (assuming |X| < 50), \Z\ < 40 from leptonic r decays, and \XZ\ < 1080 from 
B^ — !> Tu. In this work we will derive constraints on \X\ from t — >■ if ^6. 

A particularly important constraint on the mass and couplings of H^ in a 2HDM/MHDM 



is the decay 6 — )■ S7 (24, |25[, which has been measured to be in agreement with the SM 
prediction. It is the combination of couplings XY* and \Y\^ which enters the decay rate 
for b -^ s'-f (the contribution from |Xp can be neglected). Since XY* = 1 in Model II and 
the flipped 2HDM (i.e. d ^ w), the stringent bound vnii± > 295 GeV at 95% c.l can be 
derived for all values of tan/3 |26|. In contrast, in the MHDM with d ^ u, the combination 
XY* is only weakly constrained by the sum rule in eq. ([3]), and can be negative. For d = u 
(e.g. Model I and the leptonic-specific 2HDM) the constraint in eq. ([3]) does not apply. 
Therefore a light H^ (i.e. mH± < mt) being compatible with b — ^ S7 is still a possibility. 



Recent studies of the bounds on H of the MHDM from b ^ S'j [28|, |29[ derive the following 



approximate 2a intervals for the real part of XY* with mH± = 100 GeV: 

-l.KReXF* <0.7. (7) 

In deriving this constraint it is assumed that |1^| is not so big (e.g. \Y\ < 1, which is 
required from other low-energy processes such as Z — )■ bb). There is also a constraint on the 
imaginary part of XY* from a different process [IraXY* < 0.1), but for simplicity we will 
consider X and Y to be real. 




FIG. 1: The region of the [X, y, Z] plane ahowed by the unitarity constraint of eq. dH]). The 
constraint from 6 — > 57 is shown as \XY*\ < 1.1 for Re{XY*) < 0, and \XY*\ < 0.7 for Re{XY*) > 
0. 



In addition to the above phenomenological constraints, there are constraints on X, Y 
and Z from the unitarity of the matrix U for the democratic 3HDM, which were studied in 



171 . |30|. There is a non-trivial relationship on the couplings of the lightest H given as: 



A2 \U- 



Id 



+ \Y2?\Ul\ + 



Z2?\Ul 



(8) 



This constraint ensures that the magnitudes of X2, Y2 and Z2 cannot all be simultaneously 
less than one, or all be simultaneously greater than one. This is due to the fact that all 
three vacuum expectation values {vd,Vu,Ve) cannot be simultaneously large or small. 

In our numerical analysis we will always take |y| < 0.8 (as discussed above), and so 
the requirement that all three couplings cannot be simultaneously greater than unity is 
automatically satisfied. We will be concerned with the parameter space of |X| >> |y|, \Z\ 
(which corresponds to large tan/3 and moderate/small tan 7), a choice which satisfies the 
requirement that all three couplings cannot be simultaneously less than unity. In Fig. ([8]) we 
show the region in the plane [X, F, Z] allowed by the unitarity constraint in eq. ([8]), imposing 
\Y\ < 0.3 and the constraint on XY* in eq. ([7j). It can be seen that there is no parameter 
space where both of \X\ and \Z\ are less than unity, and the parameter of interest to us (i.e. 
\X\ » \Y\, \Z\) is fully compatible with the unitarity constraint in eq. ([8]). This constraint 
on the couplings is removed in a 4HDM due to the presence of a fourth vacuum expectation 
value. 



B. The Aligned 2HDM (A2HDM) 



The A2HDM is a 2HDM in which NFC is not imposed (271, and both scalar doublets ($1 
and $2) couple to all types of fermions. Tree-level FCNCs are eliminated by imposing an 
alignment of the Yukawa couplings of $1 and $2- The interaction of H^ with the fermions 
in the A2HDM can be written in the same way as eq. ([T]), but the couplings X, Y, and Z are 
determined by five parameters (instead of the four parameters in the democratic 3HDM) 
and the unitarity constraint of eq. ([8]) does not apply. Apart from these two differences 
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(which were discussed in (30|), the phenomenology of if^ in the democratic 3HDM and the 
A2HDM is essentially the same, and our numerical results will apply equally to both models. 
In particular, the magnitudes of X,Y and Z determine the BRs of if ^. In the A2HDM the 
extra free parameter can be taken to be a phase in the coupling Y, and such a phase does 
not have an effect on the BRs of H^. Moroever, we will be concerned with the parameter 
space of \X\ » \Y\, \Z\, which is compatible with the unitarity constraint in a 3HDM. 
Hereafter, when the text refers to "if ^ of the MHDM" , the implicit meaning is for an H^ 
of a 3HDM, the A2HDM, and for a MHDM with more than three scalar doublets. 

III. A LARGE BR(if ± -^ cb) IN THE MHDM AND A2HDM 

In a MHDM and in the A2HDM the expressions for the partial widths of the decay modes 
of H^ are: 

r(ii± -, tu) = G^^^-^'^^\' (9) 

47rv 2 

In r(ii^ — )■ ud) the running quark masses should be evaluated at the scale of mf{±, and 
there are QCD vertex corrections which multiply the above partial widths by (l + 17as/(37r)). 
In the 2HDM the parameter tan/3 determines the magnitude of the partial widths. The 
branching ratios are well known, and for the case of interest of mfj± < rrit one finds that 
the dominant decay mode is either H"^ — )■ cs or H^ — )■ rz/, depending on the value of tan/3. 
In model I the BRs are independent of tan/3, and BR(ii''= — )■ tu) is about twice that of 
BR(ii± -^ cs). 

The magnitude of BR(ii''= — )■ cb) is always less than a few percent in three (Models I, 
II and lepton-specific) of the four versions of the 2HDM with NFC, since the decay rate is 
suppressed by the small CKM element Vcb (^ Vcs)- In contrast, a sizeable BR(if^ — )■ cb) can 
be obtained in the flipped 2HDM for tan/3 > 3. This possibility was not stated explicitly in 
[Tj when the flipped 2HDM was discussed. The first explicit mention of a large BR(ii^ — >• 
cb) seems to have been in [l7|], and a quantitative study followed soon afterwards in |l8| . 
As discussed in section II, the condition mH± < TUt in the flipped 2HDM would require 
additional New Physics in order to avoid the constraint on mH± from b — > 57, while this is 
not the case in the MHDM. 



A. The dominance of BR(H^ -> cb) for |X| >> \Y\, \Z\ 

A distinctive signal of H^ from a MHDM for mH± < fnt would be a sizeable branching 
ratio for H^ — )■ cb. For mH± < ^t, the scenario of |X| >> |y|, |Z| in a MHDM gives 
rise to a "leptophobic" H^ with BR(ii± -^ cs)+BR{H^ -^ cb) ~ 100%. Consequently, 
BR(if='' — )• Tu) is negligible (<< 1%). The other decays of ii^ to quarks are subdominant, 
with BR(ii^ — )■ us) ~ 1% and BR(ii^ — >■ t*b) only becomes sizeable for rriHi ~ TUt, as can 
be seen in the numerical analysis in [2l| in the flipped 2HDM. Note that the case of |X| >> 
\Y\, \Z\ is obtained in the flipped 2HDM for tan/3 > 3, because \X\ = tan/3 = l/\Y\ = l/\Z\ 
in this model. 



In the scenario of \X\ » \Y\, \Z\ the ratio of the two dominant decays, BR(i7^ — > cb) 
and BR(if^ — )■ cs), approaches a constant value, which is given as follows: 



BR(f/± ^ cs) IK, 1 2m 

The CKM elements are well measured, with Vcb ~ 0.04 (a direct measurement) and Vcs ~ 
0.97 (from the assumption that the CKM matrix is unitary). The running quark masses m, 
and rrih should be evaluated at the scale Q = mH± , and this constitutes the main uncertainty 
in the ratio Ri,s- There is relatively little uncertainty for m^, with mb {Q = 100 GeV) ~ 3 
GeV. There is more uncertainty in the value of m,, although in recent years there has 
been much progress in lattice calculations of rris, and an average of six distinct unquenched 



calculations [36[ gives mg = 93.4 ±1.1 MeV [37[ at the scale oi Q = 2 GeV. A more 
conservative average of these six calculations, nis = 94 ± 3 MeV, is given in |38|. In [39(] 
the currently preferred range at Q = 2 GeV is given as 80 GeV < m, < 130 MeV. Using 
rris = 93 MeV at the scale of Q = 2 GeV (i.e. roughly the central value of the lattice 



averages in [37|, |38|) one obtains mg {Q = 100 GeV) ~ 55 MeV. Taking m, = 80 MeV and 
rris = 130 MeV at Q = 2 GeV one obtains m, ~ 48 MeV and rris ~ 78 MeV respectively at 
Q = 100 GeV. 

Smaller values of rris will give a larger BR(if''= — )■ cb), as can be seen from eq. flTTl) . Note 
that the value m, = 55 MeV is significantly smaller than the typical values rng ~ 150 — ?■ 200 
MeV which were often used in Higgs phenomenology in the past two decades. We emphasise 
that the scenario of|X| >> |y|,|Z| with mH± < rrit has a unique feature that the magnitude 
of rUg is crucial for determining the relative magnitude of the two dominant decay channels 
of H^. This is not the case for most other non- minimal Higgs sectors with H"^ that are 
commonly studied in the literature. 

In [l8| the magnitude of BR(if "*= — )■ cb) in the MHDM was studied in the plane of |X| and 
|y|, for \Z\ = and 0.5, taking rus = 0.18 GeV and rrib = 5 GeV. With these quark masses 
the maximum value is Rbs = 1.23, which corresponds to BR{H'^ -^ cb) ~ 55%. However, 
the values of m, = 0.18 GeV and m^ = 5 GeV are not realistic (as was subsequently noted in 



[20|), and two recent papers [21|, |23| have updated the magnitude of Rbs in the flipped 2HDM 
using realistic running quark masses at the scale oimjj±. In [2l|, it appears that mg = 0.080 
GeV at the scale of mjj± was used, which gives BR^H"^ -^ cb) ~ 70%, in agreement with 
our results. In [23|, rris = 0.077 GeV at the scale of mH± was used, with a maximum value 
for BR{H^ — )■ cb) of ~ 70%. We note that none of these papers used the precise average 



rris = 93.4 ±1.1 MeV [37j of the lattice calculations, which gives rris ~ 55 MeV at the scale 
of mH±- This smaller value of rris leads to a maximum value of BR{H^ — )■ cb) which is 



larger than that given in [18|, |2l|, |23| , as discussed below. 

We now study the magnitude of H^ — )■ cfe as a function of the couplings X,Y,Z. In 



Fig. ([2^) we update the numerical study of [18| for BR[H — )■ cb) in the plane [X, Y] in a 



MHDM with \Z\ = 0.1, using m, = 0.055 GeV and rrib = 2.95 GeV at the scale oimH± = 120 
GeV. With these values for the quark masses the maximum value is BR(if''= -^ cb) ~ 81% 



i.e. a significantly larger value than BR{H — )■ cb) ~ 55% in [18|. Taking a lower value 
of rris = 0.08 GeV one has BR{H^ -^ cb) ~ 69%, and for m, = 0.048 GeV one has 
BR{H'^ — > cb) ~ 86%. In Fig. ([2^) we also display the bound from b -^ sj (for mH± = 100 
GeV), which is \XY\ < 1.1 for XY* being real and negative, and \XY\ < 0.7 for XY* being 
real and positive. The parameter space for BR{H^ — >■ cb) > 60% roughly corresponds to 
|X| > 1 and |y I < 0.25 for \XY\ < 0.7. In Fig. di) and Fig. (^) we show BR{H^ -^ cs) and 
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BR(iir^ — )■ Tu) respectively. As expected, BR(if^ — ;■ cs) is maximised for \Y\ » |X|, \Z\ 
while BR(ilf^ — >■ rz/) is maximised for \Z\ » |X|,|y|. In Fig. ([3]) we show contours of 
BR{H^ -^ ch) in the plane [X, Z] for m^i = 120 GeV and |r | = 0.05. For this value of |F| 
the constraint from 6 — ;■ 57 is always satisfied for the displayed range of \X\ < 20. One can 
see that the largest values of BR(if^ -^ ch) arise for \Z\ < 2. 

B. The decay H^ — ^ A^W* for m^o < mj|^± 

The above discussion has assumed that H^ cannot decay into other scalars. We now 
briefly discuss the impact of the decay channel H"^ -^ A^W*, which has been studied in 
the 2HDM (Type II) in |40J and in other 2HDMs with small |X|, \Y\ and \Z\ in 0, with 



direct searches at LEP (assuming A^ — )■ hb) performed in 41|. In a general non-SUSY 
2IIDM the masses of the scalars can be taken as free parameters. This is in contrast to the 
MSSM in which one expects mii± ~ rriAo in most of the parameter space. The scenarios 
of rriAo < mH± and m^o > rriHi are both possible in a 2HDM, but large mass splittings 
among the scalars lead to sizeable contributions to electroweak precision observables |42| . 



which are parametrised by the S, T and U parameters [43|. The case of exact degeneracy 
{rriAo = rriHO = mni) leads to values of S, T and U which are almost identical to those 
of the SM. A recent analysis in a 2HDM [4^ sets rrijjo = ?n^o and sin(/3 — a)=l, and 
studies the maximum value of the mass splitting Am = rrijjo — mu± (for earlier studies 
see |45[). For m^io = 100 GeV the range —70 GeV < Am < 20 GeV is allowed, which 
corresponds to 80 GeV < mH± < 170 GeV. For 772^0 = 150 GeV the allowed range is 
—70 GeV < Am < 70 GeV which corresponds to 80 GeV < m//± < 220 GeV. Consequently, 
sizeable mass splittings (of either sign) of the scalars are possible. Analogous studies in a 
MHDM have been performed in [46|, with similar conclusions. 

If m^o < rnH± then the decay channel H^ -^ A^W* can compete with the above decays 
of if ^ to fermions, because the coupling H^A^W is not suppressed by any small parameter. 
In Fig. (HJi) we show contours of BR(if± -> A^W*) in the plane [X,Y] with \Z\ = 0.1, 
m^io = 80 GeV and m/^± = 120 GeV. The contours are essentially vertical in the parameter 
space of interest (i.e. \Y\ < 0.5 and \X\ » 1) because the contribution of the term m^|yp 
to the decay widths of H^ to fermions is small. Comparing Fig. (jl^) and Fig. (|2^) one can 
see that for \X\ ~ 5 both BR{H'^ — )■ A^W*) and BRi^H"^ — )■ cb) are dominant, with roughly 
equal BRs. For smaller m^o (e.g. < 80 GeV) the contour of BR{H^ — )■ A'^W*)=50% would 
move to higher values of \X\. Since the dominant decay of A^ is expected to be A^ — )■ bb, the 
detection prospects in this channel should also be promising because there would be more 
b quarks from t -)■ H^b, H^ -)■ A^W*, A° -)■ bb than from t -^ H^b with H^ —)■ cb. We 
note that there has been a search by the Tevatron for the channel t — )■ H^b, if^ — )■ yl°H^*, 
A'^ — 7- T^T^ [43], for the case of m^o < 2mi, where A'^ — )■ bb is not possible J48[. 

At present there is much speculation about an excess of events around a mass of 125 GeV 



in the search for the SM Higgs boson [3|. An interpretation of these events as originating 



from the process gg ^ A^ ^ 77 has been suggested in 49|. In Fig. (JDd) we set m^o = 125 
GeV and mH± = 150 GeV. Since the mass splitting between H"^ and A° is less than in 
Fig. (UK), the contours move to lower values of \X\, but BR{H^ -)■ A^W*)=50% is still 
possible for \X\ < 2. We note that if the excess of events at 125 GeV is attributed to a SM- 
like Higgs, then in the context of a 2HDM a candidate would be the lightest CP-even Higgs 
h^ with a coupling to vector bosons of SM strength (recent studies of this possibility can be 
found in [SOI]). This scenario would correspond to sin(/3 — a) ~ 1 in a 2HDM, and therefore 
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FIG. 2: Contours of BR(iJ=t -^ cb), BR{H^ -^ cs) and BR(if± -^ tu) in the plane [X, y] with 
\Z\ = 0.1. The constraint from 6 -^ S7 is shown as \XY*\ < 1.1 for Re{XY*) < 0, and \XY*\ < 0.7 
for Re{XY*) > 0. We take ms{Q = mH±) = 0.055 GeV and mH± = 120 GeV. 



the coupling H^h^W (with a magnitude ~ cos(/3 — a) in a 2HDM) would be close to zero. 
Hence the decay H^ — )■ h^W* would be supp ressed by this small coupling, as well as by the 
virtuality of W*. Several recent studies [SjJ, |52| fit the current data in all the Higgs search 
channels to the case of a neutral Higgs boson with arbitrary couplings. A SM-like Higgs 
boson gives a good fit to the data, although a slight preference for non-SM like couplings 
is emphasised in [52|. If the excess of events at 125 GeV turns out to be genuine and is 
well described by a non-SM like Higgs boson of a 2HDM with a value of sin(/3 — a) which 
is significantly less than unity, then BIi{H^ — )■ h'^W*) could be sizeable, with a magnitude 
given by Fig. (jU^) after scaling by cos^(/3 — a). 
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FIG. 4: Contours of BR{H^ -^ A^W*) in the plane [X, Y] with \Z\ = 0.1. The constraint from 

6 ^ S7 is shown as \XY*\ < 1.1 for Re{XY*) < 0, and \XY*\ < 0.7 for Re{XY*) > 0. In the 
left panel (a) we take mjj± = 120 GeV and m^o = 80 GeV, and in the right panel (b) we take 
mjj± = 150 GeV and m^o = 125 GeV. In both figures rUsiQ = rn^±) = 0.055 GeV. 



IV. SEARCHES FOR t 
AT THE LHC 



HH with H^ -f cs, AND PROSPECTS FOR H^ 



cb 



The case of mu± < mt + rrib with a large BR(if ^ — )■ cs) can be searched for in the decays 
of the top quark via t — )■ H^b |7|, |53[. The first discussion of t — )■ H^b followed by the decay 
H^ — )■ cb was given in [l9|. Recently, t — > H'^b with decay H^ — )■ cb has been studied in 
the context of flipped 2HDM [2l|, and in the context of the 2HDM without natural flavor 
conservation 1221]. 



There have been two dedicated searches by the Tevatron collaborations [12|, ll3| for t — )■ 
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H'^b followed by H"^ — )■ cs. The DO analysis [12| with 1 fb~^ performed a search for t — )■ H'^b 
by studying the effect of the decay H^ — )■ cs on ratios of cross sections for ti production. In 
the SM one has BR(t -)■ Wb)= 100%, and the branching ratios oi W ^ iu and W -)■ q'q 
are known. The presence of a sizeable BR(t — )■ H^b) with H^ — )■ cs would change the 
SM prediction for the ratio of the cross sections for the channels with decay W ^ iu and 
W — )■ q'q. For the optimum case of BR(if ^ — )• cs) = 100%, upper bounds on BR(t — )■ H^b) 
between 0.19 and 0.22 were obtained for 80GeV < mH± < 155 GeV. Although the decay 



H^ — )■ cs was assumed in [12| the above limits also apply (to a very good approximation) 
to the case of both H"^ — > cs and H^ — > cb having sizeable BRs, as discussed in 2l|. This 
is because the search strategy merely requires that H^ decays to quark jets. 



An alternative strategy was adopted in the CDF analysis [l3| with 2.2 fb~^. A direct 
search for the decay H^ — )■ cs was performed by looking for a peak centered at Mh± in the 
dijet invariant mass distribution, which would be distinct from the peak at M^ from the 
SM decay t — )■ Wb with W — )■ q'q. For the optimum case of BR(if ^ — )■ q'q) = 100%, upper 
bounds on BR(t — )■ H^b) between 0.32 and 0.08 were obtained for 90 GeV < mff± < 150 
GeV, with the greatest sensitivity being at mfj± = 130 GeV. No limits on BR(t — )■ H^b) were 
given for the region 70 GeV < m_f/± < 90 GeV due to the large background from W — )■ q'q 
decays. For the region 60 GeV < mH± < 70 GeV, limits on BR(t — > H^b) between 0.09 and 



0.12 were derived. As stated in [13|, the above limits also apply to other hadronic decays of 
H"^, although with slight changes in the sensitivity to BR(t — )■ H^b) because the dijet mass 
resolution depends mildly on the flavour of the quarks. The search strategy in [13|] and does 



not have sensitivity to the region 80 GeV < mH± < 90 GeV due to the large background 



from W — )■ cs, ud. The combination of the four searches at LEP for e"''e~ — ?■ H'^H~ [54 
derived the limit m^i > 81 GeV for the scenario of BR(i7''= — )■ cs) ~ 100%, with the 
following additional small intervals excluded (at 95% c.l): 86 GeV < mH± < 88 GeV and 
mH± ~ 84 GeV. Therefore the region of 80 GeV < mH± < 90 GeV with H"^ decaying 
dominantly to quarks (e.g. cs, cb) has not yet been entirely excluded yet. 

Concerning the prospects at the LHC, there has been a simulation of t — )■ H^b followed 
by the decay if^ — )■ cs by the ATLAS collaboration in [55|, assuming that one of the 



top/antitop quarks in the tt events decays leptonically via t — >■ Wb — >■ iub. This strategy 



is very similar to the CDF analysis of [13[, and directly looks for a peak centered at mH± 
in the invariant mass distribution of the jets from H^ — )■ cs. Two 6-tags are applied, and 
the peak from H^ — )■ cs is obtained by reconstructing the two untagged jets. The mass 
resolution of the peak can be further improved by full reconstruction of the ti event. For 
y/s = 7 TeV with 1 fb~^ of luminosity, values of BR(t — )• H'^b) as low as 0.04 can be probed 
for 110 GeV < mfj± < 150 GeV. This sensitivity is superior to that achieved for the decay 
t — )■ H'^b followed by if^ — )■ ru with the same integrated luminosity [l3|. Again, as in \1^ 
there is little or no sensitivity to the region 80 GeV < mH± < 90 GeV. 

The first search for t — )■ H'^b followed by the decay H^ — )■ cs at the LHC has been 



performed by the ATLAS collaboration with 0.035 fb~^ in 1^. Due to the small amount 
of integrated luminosity, only one 6-tag was applied. The limits on BR(t — )■ H^b) are 
comparable to those from the Tevatron search in [l3|, with limits of BR(t — )■ H^b) < 
0.25, 0.15 and 0.14 for mH± = 90 GeV, 110 GeV and 130 GeV respectively. 

If BR(if ^ — )■ cb) were the dominant decay channel, as can be the case in the MHDM and 



the A2IIDM, the requirement of tagging the b from H — )■ cb (as suggested in [19|, |2l|, |22| |) 
would provide sensitivity to BR(t — )■ H^b) in the problematic region 80 GeV < mH± < 90 
GeV, and should improve the sensitivity for m/f± > 90 GeV. We now estimate the gain in 
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sensitivity using realistic values for the ^-tagging efficiency (e;, = 0.5), the probability of a 
c-quark being misidentified as a 6 quark (ec = 0.1) and the probability of a light quark being 
misidentified as a 6-quark (e^ = 0.01). The two dominant backgrounds to the peak at mH± 
in the dijet invariant mass distribution are from W ^ ud and W — )■ cs, which we take to be 
equal in magnitude. For the case of BR(iir^ — )■ ch) near 80%, the ratio of the signal to the 
background {S/yB) with and without the 6-tag is given approximately as follows: 

i^^^^^~^^^^~2.13. (12) 

[5/V5]|^tag V (ej + ec) 

We encourage a detailed simulation by the Tevatron and LHC collaborations in order to 
obtain a more realistic estimate of the increase in sensitivity over the current strategy of not 
applying a 6-tag to the jets originating from if ^. 



We note that a recent paper [56[ has performed a simulation for a very similar signature 
which arises from the decay t -^ h^c — >■ bbc in a different 2HDM with FCNCs. This signature 
looks identical to the signature arising from H^ — )■ cb but there are several kinematical 
differences. The process t — )■ H^b — )■ bbc would give a peak at mi{± in the dijet invariant 
mass distribution in which only one of jets has originated from a b quark, with the other 
two 6— jets coming from the decay of ti. In contrast, for t — )■ h^c — )■ bbc both of jets in 
the dijet invariant mass distribution would originate from b quarks, while the third 6— jet 
would come from the decay of t or t. The study in |56| is specifically for t — )■ h^c — )■ bbc, 
and it was found that the sensitivity to BR(t — ?> h^c) was significantly superior to that for 
t — ;■ H^b followed by H^ -^ cs, which can be attributed to the extra 6-tag i.e. the increase 
in sensitivity in J56J compared to that obtained for the LHC simulation without the 6-tag 



m 



55[ is significantly greater than the value of 2.13 in eq. ( 1T2|) . and could be as large as a 



factor of six. 



V. NUMERICAL RESULTS 

We now quantify the magnitude of if ^ — )■ cb events produced in the decays of t quarks, 
and compare this with the expected sensitivity at the LHC. For the partial decay widths of 
t — ;■ W^b and t -^ H^b we use the leading-order expressions (with |Vtfe| = 1) as follows: 

r(t ^ WH) = ^K + 2M^][l - M^MT (13) 



T{t ^ HH) = ^^[m2|F|2 + ml\X\^][l - ml^/m: 



The multiplicative (vertex) QCD corrections to both t — > W^b and t — )■ H^b essentially 



cancel out in the ratio of partial widths 57|, and thus they do not affect BR(t — )■ H^b) 



significantly. In the phase-space function of both decays we neglect nib, and in the terms 
ml\Y\'^ and m^|Xp we use rrit = 175 GeV and rrif, evaluated at the scale of mH± (i-e. 
nib ~ 2.95 GeV). 

In Fig. ([5^) and Fig. (|5]d) we show contours of the sum of 

BR(t -> HH) X [BR(ii^ -> cs) + BR(ii^ -^ cb)] (14) 

in the plane of [X, Y] for niH± =80 GeV and niH± =120 GeV respectively, setting \Z\ = 0.1. 
The cross section in eq. (ITll) is the signature to which the current search strategy at the 
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Tevatron and the LHC is sensitive, i.e. one 6-tag (LHC [1J|) or two 6-tags (Tevatron [13|]) are 
applied to the jets originating from ti decay, but no 6-tag is apphed to the jets originating 
from H^. For the case of [BR(if^ — )■ cs)+BK{H'^ — )■ c6)]=100% the current experimental 
limits for m^i = 120 GeV are BR(t -^ H^h) < 0.14 from ATLAS with 0.035 fb^^ [13, 
BR(t -^ H^h) < 0.12 from CDF with 2.2 fb'^ [ij, and BR(t -^ H^b) < 0.22 from DO with 
1 fb~^ [12]. In Fig. (EId) for mH± =120 GeV these upper limits would exclude the parameter 
space of |X| > 40 and small \Y\ which is not excluded by the constraint from b — )■ s'y. For 



the mass region 80 GeV < mi{± < 90 GeV there is only a limit from the DO search in [12 
which gives BR(t — )■ H'^b) < 0.21. From Fig. ([5^), for mH± = 80 GeV one can see that this 
limit excludes the parameter space of |X| > 35 and small \Y\. 

In both Fig. (|5^) and Fig. ([5)d) we show contours of 1%, which might be reachable in 
the 8 TeV run of the LHC. Simulations by ATLAS (with v^ = 7 TeV) for H^ ^ cs ^ 
have shown that the LHC should be able to probe values BR(t — )■ H^b) > 0.05 with 1 fb~^ 
for rriH^ > 110 GeV, with the greatest sensitivity being around m/f± = 130 GeV. For the 
operation with ^/s = 8 TeV and an anticipated integrated luminosity of 15 fb^^ one expects 
increased sensitivity (e.g. BR(t — )■ H^b) > 0.01 for mH± > 110 GeV), although the region 
80 GeV < mH± < 90 GeV might remain difficult to probe with the strategy of reconstructing 
the jets from H^. An alternative way to probe the region 80 GeV < mH± < 90 GeV is to use 



the search strategy by DO in [12|, and presumably the LHC could improve on the Tevatron 
limit on BR(t -^ H^b) < 0.21 for 80 GeV < mH± < 90 GeV. From Fig. ([5b) (for m^i = 120 
GeV) one can see that the region of \Y\ > 0.2 and \X\ < 4, which is not excluded by 6 — )■ 57, 
would be probed if sensitivity to BR(t -^ H^b) > 0.01 were achieved. However, a large part 
of the region roughly corresponding to \Y\ < 0.2 and \X\ < 10 (which is also not excluded 
by 6 — )■ 57) would require sensitivity to BR(t — )■ H^b) < 0.01 in order to be probed with 
the current search strategy for t — )■ H^b, and this is probably unlikely in the 8 TeV run of 
the LHC. 

Increased sensitivity to the plane of [X, Y] can be achieved by requiring a 6-tag on the 
jets which originate from the decay of H^. In Figs. ([6]) and ([7]), for mfj± = 80 GeV and 
mfj± = 120 GeV respectively, we show contours of 

BR(t -^ HH) X BR(if± -^ cb) . (15) 

With the extra 6-tag, as described in eq. ( !T2|l . the sensitivity should reach BR(t — )• 
H'^b) X BR(i7^ — 7- cb) > 0.5%, and perhaps as low as 0.2%. In the latter case, one can see 
from Figs. ((Hb) and ^]p) that a large part of the regions of |X| < 5 (for mH± = 120 GeV) and 
\X\ < 3 (for mfj± = 80 GeV) could be probed, even for \Y\ < 0.2. Therefore there would 
be sensitivity to a sizeable region of the parameter space of [X, Y] which is not excluded by 
b — 7- 57, a result which is in contrast to the above case where no 6-tag is applied to the 6-jets 
originating from H^. We encourage a dedicated search for t — > H^b and H"^ — )■ cb by the 
Tevatron and LHC collaborations. Such a search would be a well-motivated extension and 



application of the searches which have already been carried out in [13| and [IJ], and would 
offer the possibility of increased sensitivity to the fermionic couplings and mass of H^ in 
the A2HDM and a MHDM. 



VI. CONCLUSIONS 

Light charged Higgs bosons (H^) are being searched for in the decays of top quarks 
{t — !■ H'^b) at the Tevatron and at the LHC. Separate searches are being carried out for 
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FIG. 5: Contours of the sum of BR(t -^ HH) x BR(i7± -^ cs) and BR(t -^ HH) x BR(if± ^ ch) 
in the plane [X, Y] with |Z| = 0.1, where mu± = 80 GeV (left panel) and mH± = 120 GeV (right 
panel) . 




FIG. 6: Contours of BR(t -^ HH)xBR{H^ -^ cb) in the plane [X, Y] with \Z\ = 0.1 for mj;^± = 80 
GeV. The constraint from 6 — )■ 57 is shown as |Xy*| < 1.1 for Iie{XY*) < 0, and |Xy*| < 0.7 for 
Re(Xy*) > 0. We take rrisiQ = '>^H±) = 0-055 GeV and show the range < |y| < 0.8 (left panel) 
and < |y| < 0.3 (right panel). 



the decay channels H^ -^ cs and H^ — ?> ru, with comparable sensitivity to the mass and 
fermionic couplings of H"^. The searches for H^ -^ cs in [13| and IJ] look for a peak at 
mH± in the dijet invariant mass distribution, with the assumption that neither of the quarks 
is a fe quark. 

In some models with two or more Higgs doublets (the Aligned 2HDM and a MHDM with 
three or more scalar doublets) the branching ratio for H^ — )■ cb can be as large as 80%. 
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FIG. 7: Contours of BR(t -^ HH) x BR(i7± -^ cb) in the plane [X, Y] with \Z\ = 0.1 for 
mj|^± = 120 GeV. The constraint from 6 — > 57 is shown as \XY*\ < 1.1 for Re{XY*) < 0, 
and \XY*\ < 0.7 for Re{XY*) > 0. We take ms{Q = mjy±) = 0.055 GeV and show the range 
< |y| < 0.8 (left panel) and < |y| < 0.3 (right panel). 



Moreover, such a H^ could be light enough to be produced via t — >■ H^b, as well as respect 
the stringent constraints from 6 — )■ S7 on both mfj± and the fermionic couplings of H"^. 
This is in contrast to H^ in other 2HDMs for which a large branching ratio for H^ — )■ cb 
is possible (such as the flipped 2HDM for mH± < rrit), but one expects mH± > frit in order 
to comply with the measured value of 6 — ?■ S7. In the context of the Aligned 2HDM and 
a MHDM we suggested that a dedicated search for t — )■ H'^b and H"^ — )■ cb would probe 
values of the fermionic couplings of if ^ which are currently not excluded by measurements 
of 6 —7- 57. Such a search would require a 6-tag of one of the jets originating from H^, 
and would afford sensitivity to a smaller value of the branching ratio of t — )■ H^b than that 
obtained in the ongoing searches, which currently do not make use of this additional 6-tag. 
We emphasised that a dedicated search for t -^ H^b and H^ — )■ cb at the Tevatron and LHC 
would be a well-motivated and (possibly) straightforward extension of the ongoing searches 
for t — )■ H^b with decay H"^ — )■ cs. 
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